Hyperbolic decoupling of tangent space and effective dimension of dissipative systems 
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We show, using covariant Lyapunov vectors, that the tangent space of spatially-extended dissi- 
pative systems is split into two hyperbolically decoupled subspaces: one comprising a finite number 
of frequently entangled "physical" modes, which carry the physically relevant information of the 
trajectory, and a residual set of strongly decaying "spurious" modes. The decoupling of the physical 
and spurious subspaces is defined by the absence of tangencies between them and found to take 
place generally; we find evidence in partial differential equations in one and two spatial dimensions 
and even in lattices of coupled maps or oscillators. We conjecture that the physical modes may 
constitute a local linear description of the inert ial manifold at any point in the global attractor. 

PACS numbers: 05.45.-a,05.70.Ln,05.90.+m,02.30.Jr 



U 



> 

in 

(N 

o 



% 



I. INTRODUCTION 

Nonlinear, dissipative, partial differential equations 
(PDEs) are frequently used for describing natural phe- 
nomena in all fields of physics and beyond. They are par- 
ticularly important for out-of-equilibrium problems such 
as pattern formation, turbulence, spatiotemporal chaos, 
etc. [ll, |2|. Because of the nonlinearities involved, their 
generic solutions are mostly studied via numerical simu- 
lation. 

A simple question arises then: can one faithfully inte- 
grate such PDEs, which are formally infinite-dimensional 
dynamical systems, when numerical schemes obviously 
involve only a finite number of degrees of freedom? For 
the practitioner, the answer is yes, as one easily ob- 
serves that, increasing numerical resolution, the obtained 
(chaotic) solutions converge, in the sense that they ex- 
hibit the same dynamical or statistical properties. Math- 
ematically speaking, the notion of inertial manifold, on 
which all trajectories evolve after a short transient, is 
crucial in this context [3|-[5| . For generic parabolic PDEs 
such as the Kuramoto-Sivashinsky (KS) equationj2|] and 
the complex Ginzburg-Landau (CGL) equation [6|, it is 
in fact proven that trajectories are first exponentially 
attracted to a finite-dimensional inertial manifold and 
eventually settle in a global attractor of finite Hausdorff 
dimension [T-Tl]. The inertial manifold is a smooth ob- 
ject embedding the global attractor, which in principle 
reduces the evolution of trajectories to a finite set of ordi- 
nary differential equations. However, this mathematical 
object remains largely formal, for lack of a constructive 
way to determine which modes actually constitute it. 

In fact, even the dimension of the minimal inertial 
manifold stays beyond the reach of mathematical argu- 
ments. So far, only upper bounds have been calculated, 
such as L^-^^ for the one-dimensional KS equation in a 



domain of size L [8|, [9| . This is to be contrasted with the 
intuitive expectation that it should grow linearly with 
L, as suggested by the extensivity of chaos }VA observed 
routinely in generic one-dimensional systems |13l-[l7| in- 
cluding the KS equation [13. 14]. 

In a different context, Cvitanovic and coworkers have 
developed a constructive approach to unravel the "skele- 
ton" of the chaotic dynamics of dissipative systems in- 
cluding PDEs |18M2Q|: by exploiting hierarchies of un- 
stable periodic orbits, chaotic solutions are represented 
rather faithfully in simplified, finite-dimensional spaces, 
and some quantitative properties can be systematically 
estimated. But it is not clear how the finiteness of the 
dimension of the inertial manifold is translated in this 
context, especially since it is well known that there are 
infinitely-many unstable periodic orbits involved, to var- 
ious degrees, in the dynamics. 

In the present paper, we show results that could con- 
stitute a first step toward a constructive approach to the 
finite-dimensional inertial manifold of dissipative PDEs. 
Using Lyapunov analysis, we study the tangent-space 
evolution of spatially-extended dynamical systems. Since 
A/'-dimensional dynamical systems have N Lyapunov ex- 
ponents, the numerical integration of a given PDF pro- 
vides as many Lyapunov exponents as one likes by just 
increasing the spatial resolution. Now, suppose that this 
PDF has an inertial manifold of finite dimension and 
that we use a resolution with a larger number of degrees 
of freedom. In a recent work [21|, we showed that co- 
variant Lyapunov vectors, which span the Oseledec sub- 
spaces [22] and give the intrinsic directions of growth of 
perturbations for each Lyapunov exponent, exhibit a de- 
composition of tangent space into a fixed, finite number 
of "physical" modes and a remaining set of "spurious" 
modes, whose number increases with increasing resolu- 
tion. The excess spurious modes are associated to very 



negative Lyapunov exponents and "hyperbolically" de- 
coupled from the physical modes. Hence, perturbations 
along them quickly decay and do not spread to any phys- 
ical mode. Below, we pursue this approach, and provide 
a precise characterization of the decoupling between the 
physical and spurious modes. We show in particular that 
the decoupling takes place even between arbitrary com- 
binations of physical or spurious modes, defining a finite- 
dimensional manifold in the tangent space, which should 
contain all the relevant information for the phase-space 
dynamics. We also provide a criterion giving an accurate 
estimate of the dimension of the physical manifold from 
finite- time simulations. 

As already mentioned above, the key quantities to ac- 
cess this hyperboHc decoupling are covariant Lyapunov 
vectors (CLVs). At each point in phase space, they con- 
stitute the intrinsic directions of growth of perturbations 
associated with each Lyapunov exponent, in other words, 
they span the Oseledec subspaces [22]. Therefore, their 
spatial structure is meaningful and their relative angles 
rule the hyperbolicity properties of the dynamical sys- 
tem. Note that they cannot be replaced, even quali- 
tatively, by the Gram- Schmidt vectors, which are by- 
products of the standard method [23, [2J ^^ compute Lya- 
punov exponents. It is only recently that CLVs have be- 
come rather easy to compute numerically, notably thanks 
to the efficient algorithm presented in \25] . In the present 
study, we compute the Lyapunov exponents by the stan- 
dard Gram-Schmidt method 1231 , |2^ and the associated 
CLVs by Ginelli's algorithm |25| . 

The paper is organized as follows. In Section [TTl we 
demonstrate our main results using the one-dimensional 
(ID) KS equation, a prototypical dissipative PDE show- 
ing spatiotemporal chaos. We then apply the same ap- 
proach to the ID CGL equation in Sec. Illll to show how 
the hyperboHc decoupling reflects phase-space dynamics 
in different regimes of spatiotemporal chaos. In particu- 
lar we study what is called the phase turbulence regime 
[6[ of the CGL equation, in order to answer the long- 
standing question regarding whether this regime can be 
fully described by "phase modes" . In Sees. IIVI and |Vl 
we study ID lattice systems and the 2D KS equation, 
respectively, demonstrating the general existence of the 
physical manifold. Section IVll contains a discussion and 
our conclusions. 



II. ID KURAMOTO-SIVASHINSKY EQUATION 
A. Definition and numerical scheme 
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We first focus on the ID KS equation pj, |2|: 
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XG[0,L], (1) 



with a real- valued field u{x^t). Boundary conditions 
are set to be periodic (PBC) u{x^ t) = u{x + L, t) for 
most of the data below, but rigid boundary conditions 
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FIG. 1. (color online). Spectrum of the Lyapunov exponents 
A*^^-* for the ID KS equation with L = 96. (a) Full spectra. 
Different symbols correspond to different cutoff wave num- 
bers kcut and different boundary conditions, either PBC or 
RBC, as indicated in the legend. Inset: Close-up around the 
threshold, (b) Close-up of (a) showing the positive end of 
the spectra, (c) The same Lyapunov spectra shown against 
j^. The data for kcut = 42 • 27t/L with PBC (black upward 
triangles) are omitted for the sake of clarity. The value of the 
Kaplan- Yorke dimension is about 21.6 for kcut — 42 • 27i/L 
and PBC. 



(RBC) 2i(0, t) = u{L, t) = are also used for compari- 
son. The size is fixed at L = 96 unless otherwise indi- 
cated, for which the ID KS equation exhibits spatiotem- 
poral chaos with both boundary conditions. For numer- 
ical integration, we use the pseudospectral method with 
discrete Fourier or sine modes up to a cutoff wave num- 
ber kcut' The number of collocation points is chosen 
so that no aliasing may occur. Integration is carried 
out with the operator-splitting method, which adopts 
the second-order Adams-Moulton method for the lin- 
ear terms and Heun's method (two-stage second-order 
Runge-Kutta method) for the nonlinear term. Most of 
the data presented in this section are obtained with PBC, 
L = 96, and kcut = 42-27r/I/, i.e., with 43 Fourier modes, 
from integration over a period of roughly 10^ after a tran- 
sient period about 10^ is discarded. 



B. Lyapunov exponents and Lyapunov vectors 

Figure [T]( a, b) shows the spectrum of the Lyapunov ex- 
ponents A^-^^ arranged in descending order for different 
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FIG. 2. Typical spatiotemporal structure of the trajectory 
u{x,t) (a) and CLVs Su^^\x,t) (b,c) in the physical and the 
spurious regions, respectively, for the ID KS equation (L = 
96, kcut = 42 • 27r/L, PBC). The index of the CLV shown here 
is j = 8 (b) and j = 46 (c). 



spatial resolutions and boundary conditions. The Lya- 
punov spectrum consists of two parts: first a smooth 
region of positive, zero, and negative exponents, and sec- 
ond a rather steep region of negative exponents arranged 
in steps of two for PBC. The two regions are separated 
by an abrupt change in slope, accompanied by the for- 
mation of a stepwise structure for PBC (inset). Remark- 
ably, the spectra for different spatial resolutions overlap 
almost perfectly, with additional exponents due to the 
improved resolution simply accumulating at the negative 
end of the second region (upward and downward trian- 
gles). Thus the threshold separating the two regions re- 
mains unchanged (here around j = 40) upon increasing 
the resolution, though its exact index cannot be defined 
unambiguously only from the spectrum and will be given 
later on a firm basis. Note also that in the second region 
the boundary conditions only change the multiplicity of 
modes: every other mode overlaps nicely for both types 
of boundary conditions. All these observations lead to a 
speculation that the Lyapunov modes in the second re- 
gion are "spurious" and dampen quickly, whereas all the 
physical properties of the dynamics are carried by the 
Lyapunov modes in the first region. This intuition shall 
be substantiated in the following by studying the CLVs 
8u^^\x^t) associated with the Lyapunov exponents A*^-^^. 
Figure [2] shows typical spatiotemporal structures of 
CLVs in the physical and the spurious regions of the Lya- 
punov spectrum as well as that of the trajectory for the 
resolution k^ut = 42 • 2t: / L with PBC (corresponding to 
the black upward triangles in Fig. [1]). At first glance, 
the CLVs for the physical and spurious Lyapunov modes 




FIG. 3. (color online). Spatial power spectrum of the 
CLVs for the ID KS equation (L = 96, /ccut = 42 • 
27r/L, PBC). (a) Power spectrum of CLVs of index j — 
1, 16, 32, 38, 44, 52, 60, 68, 76, and 84 (from left to right at the 
peak positions), (b) Peak wave number /cpeak of the power 
spectra (red circles) and k — [j /2] • 27r/L (blue line). 



have qualitatively different structures; whereas CLVs of 
the physical modes [Fig. Efb)] evolve in a manner similar 
to the trajectory itself [Fig. [Sfa)], those of the spurious 
modes are essentially sinusoidal apart from a slight mod- 
ulation of the amplitude [Fig. Efc)]. This can be seen 
more clearly in the spatial power spectra of the CLVs 
[Fig. [3l^a)]. In contrast to the smooth structure of the 
spectra found for the physical modes, those for the spu- 
rious modes are characterized by the existence of a sharp 
peak for the sinusoidal structure, whose wave number 
^peak increases linearly with the index j [Fig. [3fb)]. For 
PBC, the two modes forming a step in the Lyapunov 
spectrum are found to have the same peak wave number 
and even the same power spectrum, but with an arbi- 
trary phase shift in the vector structure. This multiplic- 
ity of two does not exist with RBC, since the phase of 
the sinusoid is fixed at the boundary. Moreover, the peak 
wave number kr^^]^ turns out to be simply the jth wave 
number, with the multiplicity taken into account, allowed 
in the given spatial geometry: kr^^]^ — [j l'^] ■ ^n/L for 
PBC [Fig. [3jb)], where [x] is the integer part of x. The 
power spectra of the spurious modes [Fig.[3fa)] show that 
they are dominated by a sinusoidal structure at this triv- 
ial wave number k^J^^^ and even more so as j is large. 
Therefore, the value of their Lyapunov exponent is deter- 
mined predominantly by the stabilizing linear term of the 
KS equation, i.e., the fourth-order derivative, and hence 
\U) ~ — (^peak)^ ^ ~J^ ^^^ large enough j, as is indeed 
confirmed in Fig. [TJc). 



C. Angles betAveen covariant Lyapunov vectors 

The sinusoidal structure of the Lyapunov vectors of 
spurious modes suggests that they are nearly orthogonal 
to each other. Unlike Gram- Schmidt vectors, CLVs allow 
us to check this directly through the angle 6>'^*'-^^ between 
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FIG. 4. (color online). Distributions Pv(^) of the angles 
between pairs of CLVs for the ID KS equation (/ccut = 42 • 
27r/L, PBC). (a) p^^O) versus 6* for L = 96. The indices of 
the pairs are indicated on the figure, (b) Pv(^) versus 1/0 
for L = 96 and for pairs (41,42), (43,44), (61,62), (81,82), 
(20,60) from right to left. The ordinate is averaged over both 
sides of the distribution. Inset: close-up of the distributions. 
(c) Angle distributions p^{0) for the pair (41,42) with varying 
system size L = 96, 97, . . . , 105 (from inside to outside in the 
figure). The distribution is obtained from shorter simulations 
for L > 99, recorded over time 10^. (d) Probability density 
Pvo of the distribution pv{0) at ^ = and tt, as a function 
of L (inset) and of 1/(L — Lc) with a given Lc (main panel). 
Pairs of the vectors are fixed at (41, 42) (black open symbols) 
and (43,44) (red full symbols). The value of Lc is Lc = 93.0 
for the pair (41,42) and Lc = 97.5 for (43,44). The circles 
and squares correspond to the recording time 10^ and 10^, 
respectively. 



vectors defined by the inner product 
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with L^-normalized vectors Jq 8u^^\x^tYdx = 1. Since 
the angle fluctuates in time, one needs to study its 
distribution function p^{0). 

The result is shown in Fig. HJ^a) for PBC. Indeed, the 
angle distributions Pv(^) between any pairs of CLVs of 
index j > 42 are peaked at 7r/2 and drop rapidly near 
and TT (green, blue, and indigo curves). This is also 
true if one of the vectors is taken in the region j < 41 
(purple curve), but when both vectors are taken from 
this region (black and red curves), the angle distribution 
spans the whole [0, tt] interval. The change in the distri- 
bution is found to be quite sharp in this example; notice 
that the indices are varied only one by one for neighbor- 
ing pairs in Fig. HJ^a) (the pairs in the same step, e.g., 
(42,43), are skipped). The transition between physical 
and spurious modes can thus be characterized by the ex- 



istence or the absence of vector tangencies, arbitrarily 
close to the phase-space trajectory. An angle distribu- 
tion bounded away from and tt marks the absence of 
such tangencies, while a finite weight at or tt indicates 
that the given pair of vectors has a finite probability of 
forming any arbitrarily small angle along the phase-space 
trajectory. Note that while two different, non-degenerate 
CLVs, will never be exactly coincident, thus tangent, at 
any point of the phase-space trajectory, they can get ar- 
bitrarily close to each other if the angle distribution has 
a finite value at or tt. In the following, we will use the 
term tangency to describe this latter occurrence. With 
this in mind, the results in Fig. |4l^a) show that the phys- 
ical modes exhibit such tangencies rather frequently. In 
contrast, the spurious modes, in addition to the negative 
values of their Lyapunov exponent, do not have any tan- 
gencies with other spurious modes (except the partner in 
the same step) nor with physical modes. In this sense, 
they are hyperbolically isolated from any other modes. 

A closer look reveals that, while at the threshold index 
the tail of the angle distribution decreases as Py{0) ~ 
exp(— const. /6>) with decreasing within our numeri- 
cal precision, showing an essential singularity at ^ = 
and TT, it decays even more rapidly for larger indices 
[Fig. HJb)]. This qualitative change of behavior pro- 
vides, in our opinion, a more accurate way of defining 
the threshold. Nevertheless, care should be taken when 
arguing about the existence or the absence of tangen- 
cies from finite- time simulations. To determine the ex- 
act threshold, one needs to perform a careful asymptotic 
analysis on the frequency of tangencies. This can be 
achieved by studying the system size dependence of the 
angle distribution py{6). Figure Sl^c) shows how Pv{0) 
changes for a given pair of CLVs, the pair (41,42) here, 
with varying system size L. The value of pv{0) at the 
tail, denoted by pvO, decreases with decreasing L, and 
from a certain system size the tail disappears [inset of 
Fig. m^d)]. A singularity analysis shows that it decays 
as pvO ^ exp[— const. /(L — Lc)], providing a threshold 
Lc = 93 for the pair (41,42) and L^ = 97.5 for (43,44) 
[main panel of Fig. Hl^d)]. This in turn allows us to de- 
termine the exact number of the threshold index, or the 
number of the physical modes, at A^ph = 43 for L = 96. 
The critical decay of the angle distribution observed for 
a fixed size L, Pv{0) ^ exp(— const. /6>), remains up to 
^ ^ 0, strictly, only at L = Lc for a given pair of indices; 
otherwise it should be eventually replaced by a faster de- 
cay or convergence to a finite value. 

Be aware that the threshold A'ph determined here does 
not coincide with the index at which the stepwise struc- 
ture of the Lyapunov spectrum is formed for PBC [j = 40 
here; see inset of Fig. HJa)]. Indeed, a closer scrutiny of 
the exponents reveals that the first few steps are actually 
slightly inclined, and hence no threshold can be unam- 
biguously defined from them. The exponents alone can 
only provide a good guess of the location of the threshold; 
instead, it must be determined from the hyperbolicity 
properties as evidenced in this section. 
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FIG. 5. (color online). Distribution Ps((/>) of the angle (j) 
between the physical and spurious subspaces for the ID KS 
equation (/ccut = 42 • 27r/L, PBC). (a) yOs(0) versus for L = 
96. The physical subspace is spanned by CLVs of index 1 to 
^ph = 43, while the rest spans the spurious subspace. Note 
that (j) is by definition equal to or smaller than 7r/2, in contrast 
to in Fig. [H Inset: the same data shown against 1/0. (b) 
Probability density pso of the distribution ps{0) at = for 
varying system size L. It is measured between two subspaces 
spanned by CLVs of indices given in the legend, and shown 
against L (inset) and 1/{L — Lc) (main panel) with the same 
Lc as in Fig. |4jd). The circles and squares correspond to the 
recording time 10^ and 10^, respectively. 



D. Angles betAveen subspaces 

In the previous section we showed the hyperbolic de- 
coupling of all the spurious modes from any physical 
mode. This, however, does not necessarily imply that 
they are also decoupled from the manifold spanned by 
the physical modes, as properly pointed out by Kuptsov 
and Kuznetsov [26], since it does not exclude the possi- 
bility of tangencies between a linear combination of phys- 
ical modes and of spurious ones. Here, following the idea 
in Ref. [26], we clarify this issue by studying the small- 
est possible angle formed by any linear combination of 
physical modes and any linear combination of spurious 
modes. It is equivalent to the minimum angle between 
the subspaces spanned by them. 

The angle between subspaces can be computed as fol- 
lows: Given the matrices Qph and Qsp defining the or- 
thogonal bases of the two subspaces, which can be eas- 
ily obtained by the QR decomposition of the arrays of 
the corresponding CLVs, the minimum angle between the 
two subspaces can be obtained simply from the singular 
value decomposition of the matrix Qp^Qsp [23, I2SI5 as 



COS( 



^i(QphQsp 



), 



(3) 



where ai is the largest singular value and the superscript 
^ denotes the transpose [29|. 

Figure [5] shows the result, which does not change the 
view presented above. The angle distribution Ps(0) is 
bounded away from zero [Fig. [5fa)], so that no tangency 
exists even between the two subspaces. The tail of the 
distribution is governed by the pairs of Lyapunov modes 
of closest indices from both groups, and hence we find 
a similar critical decay Ps(0) ^ exp(— const. /0) within 



the observed time window (inset). The system size de- 
pendence is also studied in the same manner [Fig. [5fb)], 
providing the same value of Lc as from the vector an- 
gle distribution Py{0) shown in Fig. IH^d). This implies 
that the subspace spanned by the A^ph physical modes 
is hyperbolically decoupled from the subspace of the re- 
maining spurious modes or from any other single spuri- 
ous mode. In strong contrast, inside the manifold of the 
physical Lyapunov modes, these modes are highly entan- 
gled with frequent tangencies as we have found in Fig. 
Ha). 



E. Domination of Oseledec splitting 

The absence of tangencies for the spurious modes can 
also be confirmed from another viewpoint, specifically, in 
terms of the domination of the Oseledec splitting (DOS) 
[30,^]. Roughly, DOS refers to the dynamical isolation 
of the Oseledec subspaces from each other due to the 
strict ordering of the local expansion rates. Let Xr (t) be 
the finite-time Lyapunov exponent obtained by averaging 
the local expansion rate of the jth CLV from time t to 
t + r, i.e.. 



Xi^\t)^-log\\Mt)Su^'\x,t)\l 



(4) 



where Jr{t) is the Jacobian operator for the evolution of 
the given dynamical system over time r and the CLVs 
are normalized as \\Su''^\x^t)\\'^ = J^ Su''^\x^t)'^dx = 1. 
Then, the splitting of a pair of one-dimensional sub- 
spaces, or CLVs, of indices ji and J2(> ji) is said to 
be dominated if xi^'\t) > X^J^\t) holds for all t with r 
larger than a finite tq. It has been mathematically proven 
that DOS implies the absence of tangency between the 
corresponding CLVs |30, 131|]. To quantify DOS, we de- 
fine, following [32|, 

AA(.^'^'^'2)(t) = X[^^'\t) - X[^^^\t) . (5) 

We then measure the time fraction of the DOS violation 



^ai,j2) = (e(AA^^'^'^'2)(t))), 



(6) 



where 0(z) is the step function and the brackets denote 
the time average. Note that one needs to compute the 
finite-time Lyapunov exponents Xr (t) from the CLVs, 
not from the standard method based on the QR decom- 
position [23^, z^ , because they indicate different local ex- 
pansion rates despite the same values of the long-time 
average. The right ones reflecting the physical properties 
of the tangent space, such as DOS, are those from the 
CLVs [30, 31]. 

The result is shown in Fig. [6l The DOS violation frac- 
tion Ur'^ for the pairs of the neighboring exponents 
with r = 0.2 drops sharply near the threshold and be- 
comes strictly zero for j > A^ph = 43 [Fig. [6l^a)] except 
the pairs in the same step. This split is found not only 
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FIG. 6. (color online). Violation of DOS in the KS equation 
(L = 96,A;cut = 42 • 27r/L, PBC). (a) Time fraction of the 
DOS violation, z/r '-^^ \ for pairs of the neighboring modes 
with T = 0.2. The pairs within the same step are omitted. 
Inset: i^I^^q 2 with j — 41 (black open symbols) and j — 43 
(red full symbols), shown against 1/(1/ — Lc) with the same 
Lc as in Fig. [UJd). The circles and squares correspond to 
the recording time 10^ and 10^, respectively, (b) v^^Iq^ for 



arbitrary pairs. The black color indicates v)-^ 



{31,32) 



0, i.e.. 



hyperbolically isolated pairs. The white dashed lines indicate 
the threshold A^ph = 43 between the physical and spurious 
modes. 



for the pairs of the neighbors; in fact u^^'^'^ stays zero 
for any pair that includes a spurious mode [Fig. [6l^b)], 
as expected from the fact that all the spurious modes 
are hyperbolically decoupled from the physical modes. 
The exact number of the physical modes A^ph can also 

be checked from the system size dependence of Vr'^ 
[inset of Fig.[6l^a)] in the same way as for the angle distri- 
butions of the vectors or the subspaces. All this confirms 
the absence of tangencies of the spurious modes, which 
has already been seen from the angle distributions, indi- 
cating the same number of the physical modes. 



F. What does the absence of tangencies imply? 

We now explore the implication of the absence of tan- 
gencies for the spurious modes. Suppose that we add an 
infinitesimal perturbation along an arbitrary CLV to the 
dynamics. If this is a spurious mode, the perturbation de- 
cays exponentially to zero as indicated by their negative 
Lyapunov exponents. Further, the absence of tangencies 
implies that this does not induce any perturbation along 
the directions spanned by the other Lyapunov modes. 
In contrast, perturbations along the physical modes will 
propagate to other physical modes through the tangen- 
cies between them, that is, whenever two modes get ar- 
bitrarily close to each other, and eventually induce ac- 
tivity in all the physical modes. This is true even if the 
initial perturbation was made along the physical modes 
with negative exponents, since they are directly or indi- 
rectly connected to those with positive exponents, where 
the propagated perturbation exponentially grows to af- 
fect the phase-space dynamics considerably. The situa- 
tion does not change if we initially add a perturbation 
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FIG. 7. (color online). Extensivity of the Lyapunov spec- 
trum (main panel) and effective dimensions (inset) of the 
KS equation (with PBC). Shown in the inset are the num- 
ber of the non-negative Lyapunov exponents (black circles) , 
the Kaplan-Yorke dimension (red squares) , the metric entropy 
(green diamonds, multiplied by 50), and the number of the 
physical modes A^ph (blue triangles). A^ph is estimated here 
from finite-time simulations of given sizes, so that it could 
be slightly underestimated. The dashed line indicates linear 
regression for each quantity. 



along a linear combination of physical or spurious CLVs: 
since the physical and spurious subspaces themselves are 
hyperbolically isolated (Fig.[5j), activities in the spurious 
modes decay altogether without propagating to one or a 
set of physical modes, and vice versa. In this sense, the 
dynamics corresponding to the physical modes is highly 
entangled, but completely decoupled from the decaying 
dynamics of the spurious modes. 

Moreover, we find that the number of the physical 
modes A^ph does not depend on the spatial resolution, 
provided it is high enough. This implies that even in the 
high resolution limit, or in the original PDE, the number 
of the physical modes remains the same, while infinitely 
many spurious modes can appear by increasing the res- 
olution. Therefore, the physical modes are probably as- 
sociated with the finite number of the intrinsic degrees 
of freedom describing the dynamics of the PDE. This 
leads to the conjecture that the physical modes may span 
the local linear approximation of the (minimal) inertial 
manifold. The number of the physical modes A^ph would 
then indicate the inertial manifold dimension, which is 
the smallest possible number of degrees of freedom to 
describe the dynamics of the given system. 



G. Extensivity 

The physical modes being the active degrees of freedom 
of the system, one expects that their number A^ph scales 
linearly with the system size L in systems known other- 
wise to show extensive chaos. This is indeed the case: 
in Fig. [Tithe Lyapunov spectra with their index rescaled 
by L collapse both in the physical and in the spurious 
regions with the stepwise structure retained. As sug- 
gested from this collapse, the number of physical modes 



is found to increase linearly with L, i.e., it is an extensive 
quantity, just like the other common quantities indicating 
effective dimensions of a dynamical system [22| such as 
the number of the non- negative Lyapunov exponents and 
the Kaplan- Yorke dimension (inset of Fig. [7j). Note that 
the physical manifold dimension Np^ (blue triangles) is 
the largest, in particular, larger than the Kaplan- Yorke 
dimension (red squares). This corroborates the specula- 
tion made above that the physical manifold may provide 
a tangent-space representation of the inertial manifold, 
which is a smooth manifold embedding the global attrac- 
tor. Though, clearly, some mathematical rigor is needed 
here, our result may provide an estimate of the KS iner- 
tial manifold dimension D at 



D = TVph ^ 0.43 X L, 



(7) 



which drastically lowers the existing mathematical upper 
bounds D < const, x L^-^^ [sl, [9|- Given the length scale 
27r>/2 associated with the most unstable wave number of 
the KS equation, which can be taken as the size of the 
"building blocks" of spatiotemporal chaos in this system, 
the above estimate implies that each building block con- 
tains roughly 2i^\f2 x 0.43 ~ 3.8 degrees of freedom. 

In fact, the estimate in Eq. ([7]) is obtained from the 
angle distributions and the DOS violation of finite-time 
simulations, so that the prefactor in Eq. ([7]) may be 
slightly underestimated. If we use, instead, the fact that 
A^ph increases exactly by two between two threshold sizes 
Lc = 93.0 and 97.5, we obtain TVph ^ 0.44 x L, or 4.0 
physical modes per building block of length 27r^/2. This 
estimate is remarkably close to the result of an argument 
by Sasa [33] showing that solutions of the ID KS equa- 
tion can be reconstructed, at least qualitatively, with four 
modes per building block. 



III. ID COMPLEX GINZBURG-LANDAU 
EQUATION 

A. Definition and numerical scheme 

To elucidate the generality of our finding, we now con- 
sider the CGL equation. It describes oscillatory mod- 
ulations in spatially extended continuous media at the 
most basic level, but the genericity of the phenomena 
described by the CGL equation is known to be in fact 
much larger [1, 6]. In one spatial dimension, it governs a 
complex field W{x^t) with x G [0,L] according to 



dW 



iy-(i + i/3)|iypw + (i + ia) 






(8) 



for which the phase ^(x,t) associated with the isochrone 
[2| of the local oscillation is given by 



^{x,t) =argiy(x,t) - p\ii\W{x,t)\. 



(9) 




We consider here only PBC, W{x ^ L,t) = W{x,t). 
Numerical integration is performed similarly to the 



FIG. 8. Spatiotemporal dynamics of the phase turbulence at 
a = 2.0 and /3 = —0.91 (a) and the amplitude turbulence at 
a = 2.5 and /3 = -2.0 (b) in the CGL equation with PBC. The 
system size is L = 64 for both cases. The gray scale covers 
the whole range of the values taken by |VK| and d$/dt except 
for the phase velocity d$/dt in the amplitude turbulence (b) 
showing frequent phase slips (bright and dark spots). Notice 
the different time scales between (a) and (b). 



KS equation, with the pseudospectral method and the 
operator-splitting method implemented by the second- 
order Adams-Moulton and Heun's methods. 

We study here two representative regimes of the spa- 
tiotemporal chaos of the CGL equation, namely the 
amplitude turbulence regime and the phase turbulence 
regime [o, 34, 6o\ . Both regimes correspond to the sit- 
uation where plane-wave solutions are rendered unsta- 
ble by the Benjamin- Feir instability. When varying the 
parameter values and crossing the Benjamin- Feir line 
at aj3 = — 1, one first encounters the phase turbulence 
regime; the phase field <l>(x,t) controls the spatiotempo- 
ral chaos, while the amplitude |VK(x,t)| remains rather 
quiescent, staying nonzero everywhere [Fig. [8fa)]. Be- 
yond this regime, amplitude turbulence takes place, in 
which defects associated with null amplitude and discon- 
tinuous change in the phase, in other words phase slips, 
are constantly created and annihilated [Fig. [Sl^b)]. The 
observed dynamics is chaotic both in amplitude and in 
phase for the amplitude turbulence, while in the phase 
turbulence regime, where no phase slips occur, it has usu- 
ally been considered that only the phase <I>(x,t) is active 
and the amplitude is effectively slaved to it [6] . Below we 
show, among other things, that this difference between 
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FIG. 9. (color online). Lyapunov spectrum for the CGL am- 
plitude turbulence. The stepwise structure starts at j = 71. 
The red dashed line indicates A^^^ = 1 - k'^ - 2{\W\'^) de- 
rived by Garnier and Wojcik for pure Fourier perturbations 



[36]. Inset: the same data against 
dimension in this case is about 28.6. 



The Kaplan- Yorke 



amplitude turbulence and phase turbulence is reflected 
in the structure of the physical Lyapunov modes, thus 
rooting it on a firm basis. 



B. Amplitude turbulence 

We first focus on the amplitude turbulence regime, in 
which the amplitude and the phase evolve chaotically on 
rather short time and length scales. Specifically, we use 
a = 2.5, /3 = -2.0, L = 64, and k^nt = 31 • 27r/L. Data 
are recorded over a period of about 10^ after a transient 
period of 2.5 x 10^. 

Figure [9] shows the Lyapunov spectrum in the ampli- 
tude turbulence regime. A stepwise region is present, 
as in the KS equation, but here the multiplicity of each 
step is four and the spatiotemporal evolution of CLVs 
in the stepwise region exhibits traveling waves (Fig. [TQ|) 
instead of the stationary wave found in the KS equa- 
tion [Fig. [2Ic)]. Traveling waves are the natural modes 
here, as can be readily seen from the linear stability 
analysis of the null solution W{x^ t) = 0, which gives 
exp[z(±/cx — ujt) + At] with A = 1 — k^ and uj = ak^. 
This makes the modes with the wave number k and 
—k distinguishable, albeit degenerate, and thus leads to 
the additional multiplicity of two in the CGL amplitude 
turbulence. Indeed, each CLV in the stepwise region 
shows a sharp peak in the power spectrum, located at 
the trivially determined wave number of the jth travel- 
ing wave, namely at /c = [{j + l)/4]27r/L (Fig. [TT]) . The 
values of the Lyapunov exponents therefore decrease as 



A(i) 



-j^ for large j (inset of Fig. [9]). A better 



estimate was obtained by Garnier and Wojcik [36], who 
derived A^-^^ = 1 — k^ — 2{\W\'^) for these stepwise expo- 
nents (red dashed line in Fig. [9|) assuming perturbations 
of pure Fourier modes. In passing, the degeneracy of k 
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FIG. 10. Spatiotemporal evolution of a typical CLV in the 
spurious region {j = 78) for the CGL amplitude turbulence. 
The phase component argdW^^^ is plotted. The two plots are 
from the same trajectory but at two distant periods of time, 
showing pure upward traveling waves and patches of upward 
and downward waves, respectively. 




FIG. 11. (color online). Spatial power spectrum Sj{k) of the 
CLVs for the CGL amplitude turbulence, shown as a function 
of the wave number k and the Lyapunov index j. Given that 
Sj{k) and Sj{—k) are statistically equivalent, their average 
is shown in the figure. The black dashed line separates the 
physical and spurious regions at Aph = 74. The blue solid 
line indicates k = [{j + l)/4]27r/L. 



and —k modes implies that these two modes can be in 
fact mixed up in a single vector; their CLVs are either in 
the pure k mode, in the pure —k mode, or patches of the 
two, changing their states very slowly in time (Fig. [TQ|) . 
In spite of these differences, the spurious modes in the 
amplitude turbulence are again hyperbolically isolated 
from any other modes. Figure [l2ja) shows the angle 
distribution Py{0) for pairs of the neighboring Lyapunov 
modes. Similarly to the KS equation, the functional 
form of pv{0) sharply changes at the threshold located 
at the pair (j, j + 1) with j = Aph = 74, above which 
the distributions are strictly bounded. The tail of the 
distribution near the threshold exhibits the same criti- 
cal decay pv{0) ~ exp(— const. /^) within our numerical 
precision [Fig. [I2jb)]. The absence of tangencies is also 
checked from the distribution of the angle (j) between the 
physical and spurious subspaces spanned by the vectors 
1 < j < Aph and Aph < j, respectively [Fig. [T2fc)]. The 
angle distribution Ps(0) is indeed bounded away from 
(j) = 0. Though the subspace angle (j) is smaller than the 
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FIG. 12. (color online). Distributions of the angle be- 
tween neighboring CLVs (a,b) and between the physical 
and spurious subspaces (c,d) in the CGL amplitude tur- 
bulence, (a) Distribution function p^{0) for CLV pairs 
(66, 67), (70, 71), (74, 75), (86, 87), (94, 95) from upper right to 
lower left, (b) p^{0) versus 1/0 for the data shown in (a). 
The ordinate is averaged over both sides of the distribution. 
(c) Angle distribution ps(0) between the physical subspace 
spanned by CLVs of indices 1 < j < A^ph = 74 and the spu- 
rious one formed by the remaining CLVs. (d) The same data 
shown against 1/0. 
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FIG. 13. (color online). Violation of DOS in the CGL ampli- 
tude turbulence, (a) Time fraction Ur'^ of the DOS viola- 
tion of the neighbor exponents as a function of r. The dashed 
line is a guide for the eyes showing an exponential decay, (b) 
^ui,J2j £q^ arbitrary pairs with r = 16. The black color in- 
dicates z/r ^'"^^ = 0, i.e., hyperbolically isolated pairs. The 
white dashed lines indicate the threshold Aph = 74 between 
the physical and spurious modes. 



angle for any pair of CLVs including a spurious mode, 
Ps(^) eventually decays as exp(— const. /0) for small (j) 
[Fig. [T2rd)]. confirming the absence of any tangencies. 
The spurious modes are therefore isolated from any other 
modes outside the step they belong to in the Lyapunov 
spectrum, and from any linear combination of the phys- 
ical modes, in the same way as for the KS equation. 
Concerning the DOS violation, though the time frac- 
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FIG. 14. (color online). Lyapunov spectrum A*^-'-* for the CGL 
phase turbulence. Different symbols correspond to different 
regions of the Lyapunov modes (see text). The Kaplan- Yorke 
dimension in this case is about 14.6. Inset: Lyapunov spectra 
with different spatial resolutions. 



tion Ur^'^^ measured with small r is not zero even in the 
spurious region, the same threshold A/'ph = 74 is found 
when we consider larger values of r [Fig.fTSTa)]. Contrary 

to what happens for ji, J2 < ^ph, Vt'^'^ for J2 > A^ph 
decreases faster than exponentially, though exponential 
decay may persist up to the largest r explorable near the 
threshold. This decay faster than exponential implies the 
existence of a finite r beyond which u^''^^ is zero. With 
such large r values, every quartet of the spurious modes is 
indeed hyperbolically isolated from all the other modes, 
while the physical modes including those in the first ap- 
parent step are connected to the neighbors [Fig. [TST b)]. 
All these results are consistent with what we have found 
for the KS equation and determine exactly the number 
of the physical modes, here at A^ph = 74. 



C. Phase turbulence 

We now turn our attention to the phase-turbulence 
regime [Fig. EJa)], in which the dynamics is believed to 
be governed by the phase variable [6]. We choose here 
a = 2.0, P = —0.91, and L = 96. The cutoff wave 
number is kcut = 31 •27r/L unless otherwise indicated. In 
this phase-turbulence regime, the global phase difference 
is conserved because of the absence of phase slips. Here 
we set it to be zero, in which case phase dynamics is 
well-described by the KS equation, at least close to the 
Benjamin- Feir line [6]. Data are recorded over a period of 
about 10^ after a transient period of 5 x 10^ is discarded. 

In this regime, the Lyapunov spectrum shown in Fig. 
[T4l reveals an interesting structure; now there are much 
less physical modes (black circles) than in the amplitude 
turbulence and they are followed by three distinct regions 
of spurious modes. The modes in the first spurious region 
resemble those in the KS equation (red squares), while in 
the last part of the spectrum we find a structure similar 
to that in the CGL amplitude turbulence (blue triangles). 
In between, the spurious modes appear to be mixed up 
with the modes that would take part in the physical re- 
gion in the amplitude turbulence regime, as suggested 
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FIG. 15. (color online). Spatial power spectrum Sj{k) of the 
CLVs for the CGL phase turbulence, shown as a function of 
the wave number k and the Lyapunov index j. Given that 
Sj{k) and Sj{—k) are statistically equivalent, their average is 
shown in the figure. The peak wave number for each j is indi- 
cated by the black rectangles. The black dashed lines separate 
the four regions of the Lyapunov modes. The blue solid lines 
indicate k = [j/2]27r/L (upper line) and k = [{j + l)/4]27r/L 
(lower line), which are the trivial peak wave numbers found 
for the spurious modes in the KS equation and in the CGL 
amplitude turbulence, respectively. 



from the rather smooth structure of the spectrum (green 
diamonds). We shall therefore refer to these spurious 
modes as KS-like, amplitude-turbulence-like (AT-like), 
and mixed, respectively, as will be justified later. The ex- 
act indices for all the thresholds will be given as well from 
their hyperbolicity properties. Increasing the spatial res- 
olution merely results in adding further AT-like spurious 
modes at the end of the spectrum without changing the 
essential features of the existing modes (inset of Fig.[T4|). 

Figure [15] shows the spatial power spectra of the asso- 
ciated CLVs. This justifies the above classification of the 
spurious modes; in the physical region the power spec- 
trum of the CLVs reflects that of the trajectory (not 
shown) and thus it hardly changes with varying index, 
except for those associated with null Lyapunov exponents 
or too close to the threshold. For the KS-like and AT- 
like spurious modes, the CLVs are nearly sinusoidal, as 
indicated by a sharp peak in their power spectrum. Their 
peak wave number /cpeak is determined by the same trivial 
geometrical rules as those for the KS equation and for the 
amplitude turbulence, respectively: /cpeak = [j/2]27r/L 
for the former and /Cpeak = [(j + l)/4]27r/L for the latter 
(blue solid lines in Fig. [T5]) . In contrast, the mixed spu- 
rious region consists of two branches; one continues from 
the KS-like spurious modes, corresponding to the upper 
branch in the peak wave number, and the other contin- 
ues from the physical modes, forming the lower branch. 
The two branches finally merge, marking the start of the 
AT-like spurious region. 

The correspondence of these modes to the spurious 
modes of the KS equation and the CGL amplitude turbu- 
lence is not found only in the peak wave number. Figure 
[16] displays the spatiotemporal structure of a typical KS- 
like and AT-like spurious mode, showing specifically the 
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FIG. 16. Spatiotemporal evolution of a typical KS-like and 
AT-like mode in the CGL phase turbulence, (a) Amplitude 
\W\ of the dynamics. (b,c) Phase component 5$*^-^'* of the 
CLV for the KS-like spurious mode at j = 28 (b) and for the 
AT-like one at j = 91 (c). 



local phase shift d^^^'' onto the phase field ^ [Eq. (J9j)] 
induced by the corresponding CLV. It is defined by the 
relation 

SW^^^ = [SA^^^ + i{AS^^^^ + /3M(^'))]e^(^+^^"^), (10) 

obtained with A{x,t) = \W{x,t)\ and W = Ae^(*^+^inA)^ 
and is a natural counterpart of the vector component 
Su'^^^ for the KS equation. In this phase shift field, the 
KS-like spurious modes [Fig. [TCT b)] reveal essentially sta- 
tionary wave structure with amplitude modulated by the 
dynamics, in the same way as those for the KS equa- 
tion [Fig.[2]^c)]. Similarly, the same traveling wave struc- 
ture as in the amplitude turbulence regime (bottom in- 
sets of Fig. [9]) is found for the AT-like spurious modes 
[Fig.[l6]:c)]. 

In spite of this rather complicated internal structure 
of the spurious region, the essential feature of the spuri- 
ous modes remains the same; they are all hyperbolically 
isolated from the physical modes. Figure [TTJ^a-d) shows 
the angle distributions for pairs of neighboring Lyapunov 
modes near the thresholds. For the threshold between the 
physical and KS-like spurious modes, the first two visible 
doublets in the Lyapunov spectrum (24 < j < 27) belong 
in fact to the physical region [see the pair (25, 26) in Fig. 
[iTt^a.b)]. whereas the first KS-like spurious mode is found 
at j = 28 [pair (27,28) in (a,b)]. The KS-like spurious 
modes are hyperbolic with respect to any other modes 
outside the step including other KS-like doublets [pair 
(29, 30) in (a,b)] and the mixed modes [pair (33, 34) in 
(a,b)]. In contrast, the mixed spurious modes, despite the 



11 




(b) 



^Pv(e) 


1 1 1 1 

(25,26) D (29,30) 
o (27,28) ' (33,34) 
A (34,35) _ 


o 


, , , 1/e 



71/4 71/2 371/4 71 

(c) 



1.5 
(d) 



2.5 



10" 



10"^ 



10 ■ 



Pv(e) 


(90,91) 


, (85,86) J 


HN 


(86,87) 

1 


ImK' 



(0) o (82,83) (86,87) 

v^ ^ D (85,86) A (90,91) 



71/4 7l/2 371/4 71 

(e) 




10^- 



10^ 



10'^ 



10" 



Pg(.9/ PvsK+M+A 



P+K vs M+A 



PvsM P+K+MvsA" 
M vs A 




10 



10 



10"^ 



/^\ o P vs K M vs A 

PoVY/vPvsM aP+K+MvsA 



□ P vs K+M+A <l P+K vs M+A- 



nil'" 




FIG. 17. (color online). Distributions of the angle between 
neighboring CLVs (a-d) and between pairs of subspaces (e,f) 
in the CGL phase turbulence, (a) Distribution function pv{0) 
for CLV pairs (34, 35), (25, 26), (33, 34), (27, 28), (29, 30) from 
outermost to innermost. They determine the thresholds for 
the physical, KS-like, and mixed regions, (b) pv{0) versus 
1/0 for the data shown in (a). The ordinate is averaged 
over both sides of the distribution, (c) pv(0) for CLV pairs 
(82, 83), (85, 86), (86, 87), (90, 91) as indicated by the labels. 
They are relevant to the threshold between the mixed and 
AT- like regions, (d) pv{0) versus 1/0 for the data shown in 
(c). The ordinate is averaged over both sides of the distribu- 
tion, (e) Subspace angle distributions Ps(0). The abbrevia- 
tions P, K, M, and A refer to the physical modes (1 < j < 27), 
the KS-like spurious modes (28 < j < 33), the mixed spuri- 
ous modes (34 < j < 86), and the AT- like spurious modes 
(87 < j < 126), respectively. The rightmost two curves over- 
lap almost perfectly, (f) The same data as (e) shown against 
1/0. 



absence of tangencies with any physical, KS-like, and AT- 
like modes, do have tangencies inside the mixed region 
[pair (34, 35) in (a,b) and pairs (82, 83), (85, 86) in (c,d)]. 
The asymmetry (with respect to 7r/2) of some of the dis- 
tributions is due to the finite length of the simulations, in- 
dicating that correlations relax very slowly for these pairs 
of CLVs. In the AT-like spurious region, every quartet 
is again isolated from any other modes including other 
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FIG. 18. (color online). Violation of DOS in the CGL phase 
turbulence, (a) Time fraction z^r '"^ of the DOS violation for 
pairs of neighboring Lyapunov exponents as a function of r. 
Pairs with and without tangencies in the angle distributions 
[Fig. [TTT a.b)] are denoted by open and solid symbols, respec- 
tively. All pairs in the KS-like spurious region (28 < j < 33 
excluding those in the same doublet) perfectly satisfy DOS, 
i.e., i/^-^^'^'^^ = 0, even with such a small value of r as 0.2. 
(b) i/!f^^'-^^^ for arbitrary pairs with r = 4. The black color 



indicates ui- 



{31,32) 



0, i.e., hyperbolically isolated pairs. The 



white dashed lines indicate the thresholds of the four regions. 



quartets in the AT-like region [pairs (86, 87), (90, 91) in 
(c,d)]. From all these observations, we precisely deter- 
mine the number of the physical modes at A^ph = 27, 
as well as the index of all the internal thresholds of the 
spurious region as illustrated in Fig. [lH 

The violation of DOS corroborates this view (Fig. [18]). 
Though the time fraction of the DOS violation Vr'^'^'^^ is 
found to decrease sometimes not monotonically with re- 
spect to r for mixed spurious modes, presumably related 
to their long-lasting correlation, Ur^'-^'^^ decays faster 
than exponentially for any pair which lacks tangencies 
in the angle distributions [Fig. [TST a)]. The matrix rep- 
resentation of i^r ^'"^^ with a sufficiently large r, which is 
only r = 4 here, clearly shows that all KS-Hke doublets 
and AT-like quartets are completely isolated, while the 
mixed spurious modes are densely connected inside the 
mixed region but with none of the modes outside [Fig. 

mb)]. 

Similarly to the cases studied above, not only pairs 
of single Lyapunov modes, but any linear combination 
of them is strictly hyperbolic if they are chosen from 
different regions. This is best demonstrated by the an- 
gle distributions for subspaces spanning the four regions 
[Fig. fTTT eJ)]. Subspaces formed by the Lyapunov modes 
spanning different regions have always nonzero angles be- 
tween each other. The angle distribution Ps(^) shows 
the critical decay with essential singularity, i.e., Ps(^) ~ 
exp(— const. /(/)), or decays even faster [Fig. fTTTf)]. In 
particular, all the spurious modes, including mixed ones, 
do not have any tangency with any of the physical modes 
and their Lyapunov exponents are always negative. The 
two essential properties of the spurious modes are thus 
satisfied by all of them. This implies that they exponen- 
tially decay to zero without interacting with any phys- 
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ical modes, though some spurious modes are internahy 
entangled, and thus they are not expected to affect the 
dynamics in the global attractor. 

These results for the amplitude turbulence and the 
phase turbulence lead us to interpret that amplitude 
degrees of freedom are slaved in the phase turbulence 
regime, and this produces the mixed spurious region in 
place, being mixed with the KS-like spurious modes. To 
test this speculation, we compare the structure of the 
CLVs with that expected under the assumption of the 
amplitude slaved by the phase. The dynamics in the 
phase turbulence regime indicates that the amplitude 
A{x^t) = \W{x^t)\ is almost uniquely determined by the 
local phase curvature ^(x,t) = d'^^/dx^^ i.e.. 



(a) trajectory 



^ (b) vector (P), / = 8 

.10^ 0.1 1 ' 



A(x,t)^A(^(x,t)), 



(11) 



as shown in Fig. HW aV This then impHes the following 
amplitude-phase relation for the CLVs, as long as they 
result from the slaved amplitude: 



dA. 



5A^^Hx,t)^—5^^=\x,t), 



(12) 



where 5A and 5^ are the vector components corre- 
sponding to A and ^, obtained through Eq. (p!Q|) and 
(5^0) = a^(5$(-^')/Sx^, and the subscript "p" denotes the 
estimate under the assumption (pT]) . Comparing (5Ap 
with the actual amplitude component 5A^^^ , one can ex- 
amine to what extent the assumption of the slaved am- 
plitude holds for each Lyapunov mode. 

The result is shown in Fig. [TOTb-f). which displays the 

density plots of 5A^^ versus 5A^^^ for a typical CLV in 
each region of the Lyapunov spectrum. For the CLVs 
in the physical region [Fig. [TOT b)] the prediction of Eq. 
([12]) holds as accurately as the slaved amplitude assump- 
tion (fTTj) does for the dynamics [Fig. [TOTa)]. It is also 
the case, albeit with somewhat larger deviations, for the 
KS-like spurious modes [Fig. fTOT c)]. This is quite natural 
because the physical modes should reflect the dynamics, 
which is governed by the phase here, and the KS-like 
spurious modes are present in the KS equation, which 
phenomenologically describes such phase instability. On 
the other hand, the mixed spurious modes, which com- 
prise degrees of freedom not obeying the slaved ampli- 
tude assumption, exhibit very different structures. For 
those which belong to the lower branch in the peak wave 
number spectrum [Fig. [TOTd)]. the amplitude 5 A varies 
widely, while the phase component stays essentially qui- 
escent, so does 5A^. This implies that these modes com- 
prise mostly the amplitude degrees of freedom, which un- 
derpins the view that independent amplitude degrees of 
freedom are lost and become spurious in this phase tur- 
bulence regime. In contrast, the mixed spurious modes 
in the upper branch are still influenced by Eq. ([12]) [Fig. 
[TOTe)]. This indicates that these modes essentially orig- 
inate from the KS-like modes, but the slaved amplitude 
assumption holds only very loosely, probably because of 
the mixture with the spurious amplitude degrees of free- 
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FIG. 19. (color online). Density plots showing the structure 
of the trajectory (a) and of the CLVs (b-f) in the CGL phase 
turbulence, (a) Density plot for the amplitude A and the 
phase curvature ^ of the trajectory. The white dashed curve 
indicates the averaged value of A for each bin of ^ . (b-f) 
Density plots for the CLVs, comparing the actual amplitude 
component 5 A with that expected from the assumption of the 
slaved amplitude, dA^ [Eq. (|12p ]. The white dashed lines in- 
dicate 5Ap — 5 A. The abbreviations P, K, M, and A denote a 
physical, KS-like, mixed, and AT-like mode, respectively. The 
upward and downward arrows next to M denote the upper 
and lower branches, respectively, in the plot of the peak wave 
number [Fig. dU^b)]. All the data shown here are recorded 
during two periods of length 10^ separated by an interval of 
roughly 10^ 



dom in the lower branch. Finally, for the AT-like spu- 
rious modes the slaved amplitude assumption does not 
work any more [Fig. [TOTf)]. The two amplitudes rather 
rotate in the 5A-3A^ plane irrespective of the dynamics, 
which is in line with their trivial traveling wave structure 

[Fig.nac)]. 

In summary, we have found that in the phase turbu- 
lence regime the tangent-space dynamics consists of the 
physical modes and the three sets of spurious modes, all 
of which are hyperbolically decoupled from any physical 
mode. These physical modes essentially stem from the 
phase degrees of freedom of the CGL equation, unlike in 
the amplitude turbulence where both phase and ampli- 
tude participate in the physical modes. The loss of the 
amplitude degrees of freedom considerably reduces the 
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number of the physical modes and produces instead the 
mixed spurious modes. This view is supported by the 
fact that the mixed spurious region starts at the expo- 
nent value close to —2 [Fig. [MT a)]. which is the stability 
of the amplitude of a single Ginzburg-Landau oscillator. 
Our results, therefore, corroborate the qualitative pic- 
ture of the phase turbulence that the dynamics is gov- 
erned by the phase, whereas, practically, the amplitude 
modes do not serve as independent degrees of freedom 
as they are essentially slaved by the phase. It should be 
noted, however, that the physical and spurious modes do 
not correspond exactly to the phase and amplitude de- 
grees of freedom, respectively, as can be seen from the 
weak violation of Eqs. (pTj) and ([12]) in the dynamics and 
the physical modes [Fig. fToTa.b)]. This implies that the 
phase reduction is not exact in the phase turbulence, in 
line with results from past studies [S^- To extract a 
proper set of global variables describing the dynamics is 
a difficult problem, which will be briefly discussed in Sec. 

ED 



IV. ID LATTICE SYSTEMS 

So far we have studied two generic dissipative PDEs 
and found the separation of a finite number of physical 
modes from the remaining spurious modes. We have ar- 
gued that the physical modes correspond to the active 
degrees of freedom necessary to describe the dynamics 
and that they probably provide a tangent-space repre- 
sentation of the inertial manifold. Although the concept 
of the inertial manifold is of great importance in such 
PDE systems, the existence of such a smooth invariant 
manifold of dimension smaller than the phase space is 
not restricted to PDEs. We therefore study spatially dis- 
crete systems in this section, one defined with discrete 
time (i.e., map) and another with continuous time (i.e., 
ordinary differential equation). 



A. ID Lattice of Coupled Tent Maps 

First we consider the following one-dimensional lattice 
of diffusively-coupled tent maps: 

xt^' = f{xD + ^ifixU,) - 2f{xl) + f{xU)] 



i.f. 



m) 



(13) 



with i = 1,2,- •• ,L, f{x) = 1 — /i|x|, and the discrete 
Laplacian operator 

Vfixl) = fixU,) - 2fixl) + f{xU). (14) 

In such systems, obviously, varying the spatial resolution 
does not make sense. Instead, we increase the coupling 
constant K so that local dynamical variables x\ possess 
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FIG. 20. (color online). Lyapunov spectrum for the cou- 
pled tent maps ([T3|) with varying coupling constant K. Inset: 
close-up of the end of the spectrum for K = 0.65. The green 
lines show the spectrum estimated under the Fourier mode 
approximation (see text). The Kaplan- Yorke dimension for 
K = 0.65 is about 39.9. 
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FIG. 21. Spatiotemporal evolution of the dynamics xj (a), of 
a physical Lyapunov mode at j = 1 (b), and of a spurious 
Lyapunov mode at j = 251 (c) in the coupled tent maps 
with K = 0.65. The spatial profiles are plotted every 8 time 
steps in order to capture the internal structure of the period-4 
collective behavior. 



stronger correlations with their neighbors, possibly re- 
ducing the number of effective degrees of freedom. 

Here we use fi = 1.1, L = 256, and PBC. With 
this value of /i, the local map exhibits four-band chaos. 
We start from random initial conditions distributed only 
within one of the bands, and hence the coupled-map 
lattice shows period-4 collective behavior. Data are 
recorded over roughly 5 x 10^ time steps after a tran- 
sient of nearly the same length. 

Figure [2Ql shows the Lyapunov spectra of this coupled- 
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FIG. 22. Spatial power spectrum of the CLVs for the tent 
coupled-map lattice with K = 0.65, shown as a function of 
the (integer) wave number k and the Lyapunov index j. The 
dashed line indicates k = 87, which gives the smallest eigen- 
value for the couling operator 1 + [K/2)T> (see text). The 
panel (b) is a close-up of (a) near the peak wave numbers of 
the CLVs at the end of the Lyapunov spectrum. 



map lattice for two values of the coupling constant K. As 
expected, while for moderate strengths of the coupling 
the Lyapunov spectrum does not show any splitting of 
the structure (e.g., red squares in Fig. [20]), strengthening 
the coupling induces a stepwise structure at the end of the 
Lyapunov spectrum [K = 0.65 in Fig. [20] (black circles)], 
similarly to what happens in PDEs. 

Concerning the structure of the CLVs, the Lyapunov 
modes in the smooth region, particularly at the begin- 
ning of the spectrum, tend to be localized at one of the 
antinodes (or a few of them) of the wavy structure of 
the dynamics [Fig. [2lTa,b)]. Over long time scales the 
weight shifts from an antinode to another (or from some 
to others). In constrast, those in the stepwise region ex- 
hibit stationary wave structures [Fig.[2Trc)] similar to the 
spurious modes in the KS equation [Fig. [2]^c)]. 

A difference from the PDE systems is found however 
in the way the peak wave number of the spurious modes 
is determined. Figure [22] shows the spatial power spec- 
tra of the CLVs. While in PDEs the peak wave num- 
ber in the stepwise region increases monotonously ac- 
cording to trivial geometrical rules (Figs. [3] and [TS] ex- 
cept for the special case of the mixed spurious modes 
in the CGL phase turbulence), in the coupled-map lat- 
tice it forms two branches, one increasing and the other 
decreasing, which eventually join at j = L at a cer- 
tain wave number [dark red spots in Fig. [22r b)]. In 
fact, this is again a simple consequence of the nearly si- 
nusoidal structure of the CLVs in the stepwise region. 
The Fourier modes 6Wi = sm{27rik/L) and cos{27rik/L) 
with integer k are eigenvectors of the coupling opera- 
tor 1 + {K/2)V in Eq. ([13]). Their associated eigen- 
value is 1 — K[l — cos(27r/c/L)], whose contribution to 
the Lyapunov exponent is log|l — K[l — cos{27rk/L)]\. 
Since it is not a monotonous function of /c, nor is the 
peak wave number in the power spectra. For K = 0.65, 
k = 87 gives the eigenvalue closest to zero, and there- 
fore is the peak wave number for the CLV of the small- 
est Lyapunov exponent (dashed line in Fig. [22]) . More- 
over, approximating the contribution from the local map 
f{x) to the Lyapunov exponent by log/i, we obtain 
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FIG. 23. (color online). Hyperbolicity properties for the 
tent coupled-map lattice with K = 0.65. (a) Time frac- 
tion z/r ^'-^^^ of the DOS violation for arbitrary pairs with 
r = 2. The black color indicates Ur^'^"^ = 0, i.e., hy- 
perbolically isolated pairs. The white dashed lines indicate 
the threshold A^ph = 250 between the physical and spurious 
modes, (b) Close-up of (a) near the end of the spectrum. 
(c) Angle distributions pv{0) for neighboring CLVs of indices 
(246, 247), (248, 249), • • • , (254, 255) from upper left to lower 
right, (d) pv{0) against 1/0. The ordinate is averaged over 
both sides of the distribution, (e) Angle distribution Ps(0) be- 
tween the physical subspace (1 < j < 250) and the spurious 
one (251 <j< 256). (f) ps(0) against 1/0. 



A = log|l — K[l — cos(27r/c/L)]| + log/i, which indeed 
captures the qualitative structure of the observed Lya- 
punov spectrum (Fig. [20]). 

Despite this apparent difference from PDEs, we iden- 
tify the same tangent-space decoupHng into the physical 
and spurious modes in this coupled-map lattice. Mea- 
suring the time fraction i/^''^^ of the DOS violation, 
we find that the last three doublets of Lyapunov modes 
are hyperbolically isolated from all the other modes [Fig. 
[23lfa,b)]. The splitting into the physical and spurious 
modes can also be confirmed from the angle distribu- 
tions between CLVs [Fig.[23l^c,d)], as well as that between 
the physical and spurious subspaces defined by the cor- 
responding CLVs [Fig. [23l^e,f)]. If the manifold of the 
physical modes is a local approximation of the inertial 
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manifold as we speculate, our results imply that dissipa- 
tive coupled map lattices may also have an inertial man- 
ifold, whose dimension is lower than that of the phase 
space. 



B. ID Lattice of Ginzburg-Landau Oscillators 

Now we study another example of the spatially dis- 
crete dynamical systems, here defined with continous 
time. Specifically, we take the system studied recently 
by Kuptsov and Parlitz [38|, namely a one-dimensional 
lattice of CGL oscillators, 



dWj 

"of 



W^-{l^if3)\W^fW,^{l^ia) 






(15) 



with complex variables Wi(t) and i = 1,2,--- ,L. This 
can be regarded as a spatial discretization of the CGL 
equation (|8]), where the parameter h denotes the lattice 
constant, and thus behaves similarly unless h is too large. 
For small /i, the system is a good approximation of the 
CGL equation, and hence the decoupling of the physical 
and spurious modes takes place in the same manner as 
we have studied in Sec. Hill Increasing /i, besides degrad- 
ing the correspondence to the CGL equation, amounts 
to increasing the chain length hL and thus leads to an 
increase in the number of physical modes. Moreover, 
Kuptsov and Parlitz [38] reported that, as the number 
of physical modes increases with h at the positive end of 
the Lyapunov spectrum, another set of modes appears 
at the negative end, which are hyperbohcally connected 
within themselves and whose number also increases with 
h. Isolated steps of the spurious modes being replaced 
from both ends, the authors argued that the two groups 
finally meet each other and the strict hyperboHc decou- 
pling between them is lost for large h. However, their 
claim is based solely on the observation of the structure 
of the CLVs and the time fraction of the DOS violation 
with a fixed r [32., ,39\ , which cannot fully determine hy- 
perbolic decoupling. Therefore, here we revisit this prob- 
lem, performing a complete analysis to clarify the nature 
of the hyperbolic decoupling in this system, and explains 
this rather peculiar behavior reported by Kuptsov and 
Parlitz. 

In the following, we take the same parameter values as 
Kuptsov and Parlitz [38|, namely a = —2.0 and /3 = 3.0 
which correspond to the amplitude turbulence regime (cf. 
Sec. IIIIBp . and vary the value of h. We fix the system 
size at L = 32 and adopt PEG, VF^+l = Wi. Kuptsov 
and Parlitz used instead a no-fiux boundary condition 
in Ref. [38], but this difference in the boundary condi- 
tions intervenes practically only in the multiplicity of the 
stepwise structure of the Lyapunov spectrum, as we have 
seen in Sec. [Ill Time integration is performed with the 
fourth order Runge-Kutta method, over a period of about 
2 X 10^ after a transient of 2.5 x 10^ is discarded. 

Figure [24l shows the Lyapunov spectra of the CGL lat- 
tice for varying h. As reported in Ref. [381], for small h the 




FIG. 24. (color online). Lyapunov spectrum for the CGL 
lattice ([15]) with varying h. The numerically obtained spectra 
are shown by the symbols, while the black line indicates the 
spectrum of the contribution from the coupling term under 
the Fourier mode approximation for h = 0.5 (see text). The 
Kaplan- Yorke dimension is about 8.8, 11.6, and 14.7 for h = 
0.5, 0.65, and /i = 0.8, respectively. 




FIG. 25. (color online). Spatial power spectrum Sj{k) of 
the CLVs for the CGL lattice with h = 0.5 (a) and h — 
0.8 (b), shown as a function of the (integer) wave number k 
and the Lyapunov index j. Given that Sj{k) and Sj{—k) are 
statistically equivalent, their average is shown in the figure. 
Note that the color plots are in different linear scales for the 
sake of clearity. 



spectrum exhibits the stepwise structure (black circles), 
but it arises in the middle of the spectrum, as opposed to 
the CGL equation for which the steps continue up to the 
end of the spectrum. In the CGL lattice, they are instead 
sandwiched by two smooth regions. With increasing /i, 
both of these smooth regions grow (green diamonds in 
Fig.[24|) and finally merge (red squares), replacing all the 
steps in the spectrum. The power spectra of the CLVs in 
Fig. [25l show that the steps here are again due to the si- 
nusoidal structure of the CLVs, as evidenced by the sharp 
peaks for the stepwise modes in Fig.[25fa), while CLVs in 
the two smooth regions do not have obvious wavy struc- 
ture. 

Analyzing hyperbolicity properties of this system re- 
veals that, for small /i, the Lyapunov spectrum splits into 
three regions, basically according to its stepwise structure 
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FIG. 26. (color online). Hyperbolicity properties for the CGL 
lattice with h — 0.5. (a) Angle distributions Pv(^) for neigh- 
boring CLVs of indices (22, 23), (46, 47), (50, 51), and (48, 49) 
from innermost to outermost curve, (b) p^(&) against XjQ. 
The ordinate is averaged over both sides of the distribution. 
The order of the curves is the same as in (a) from lower 
left to upper right, (c) Angle distribution Ps(0) between 
the physical subspace (1 < j < A^ph) and the spurious one 
(ATph+l < j) with ATph = 22 (black solid line). The red dashed 
line shows the angle distribution between the subspaces of in- 
dices 1 < i < 46 and 47 < j < 64, showing the hyperbolic 
isolation of the tangled spurious subspace. (d) The same plots 
as (c) but shown against 1/0 (black circles for the physical and 
spurious subspaces, and red squares for the other), (e) Time 
fraction Ur'^ of the DOS violation for pairs of neighbor- 
ing Lyapunov exponents as a function of r. (f) Time fraction 
^ui,j2) q£ ^Yie DOS violation for arbitrary pairs with r = 3. 
The black color indicates z^r ^'"^^ = 0, i.e., hyperbolically iso- 
lated pairs. The white dashed lines indicate the threshold 
A^ph = 22 between the physical and spurious modes, and that 
between the physical and stepwise spurious modes 1 < j < 46 
and the tangled spurious modes 47 < j < 64. 



(Fig. [26]) . The first gap in the stepwise region already 
exhibits clear hyperbolic decoupling of the two modes at 
the edge, namely those of indices 22 and 23; for them the 
CLV angle distribution Py{0) is bounded with a decay 
faster than exp(— const. /^) in the tail [black, innermost 
line and symbols in Fig. [26l^a,b)], and the time fraction 
of the DOS violation Ur ' decreases faster than ex- 
ponentially [black circles in Fig. [26r e)]. This fixes the 




FIG. 27. (color online). Hyperbolicity properties for the CGL 
lattice with h = 0.8. (a) Angle distributions pv{0) for all 
pairs of neighboring CLVs. Indicated by the arrows are the 
distributions for the pairs including either or both of the two 
Lyapunov modes with A^-^-* = 0. All the distributions except 
the one for the two zero Lyapunov modes, shown by the bold 
arrow, have finite densities at ^ = and tt. (b) Time fraction 
^ui,j2) q£ ^Yie DOS violation for arbitrary pairs with r = 
5. The black color indicates Ur'^^ = 0, i.e., hyperbolically 
isolated pairs. 



number of physical modes at A^ph = 22. The decou- 
pling holds similarly for any pair of CLVs sandwiching 
this threshold [Fig. IMT f)], as well as between the physi- 
cal and spurious subspaces, defined at this index [black 
solid line and circles in Fig. IMTcd)]. 

The spurious subspace in this system, however, con- 
sists of the central stepwise region and the smooth region 
at the end of the Lyapunov spectrum (Fig. [M|) . In the 
former, quartets of Lyapunov modes are isolated from all 
the other modes [Fig. [2611)] as is usually the case for the 
spurious modes, while in the latter, specifically here for 
indices larger than 46, the spurious modes are internally 
connected [Fig. [26l^a,b,e,f )] in a way similar to the mixed 
spurious modes in the CGL phase turbulence [see, e.g., 
Fig. [TSTb)]. These are the modes found by Kuptsov and 
Parlitz [38], though in their work there remained ambi- 
guity in the hyperbolicity properties as discussed above. 
Here we confirm that these modes are indeed connected 
internally, but never with any physical or stepwise spuri- 
ous mode, both from the DOS assay with varying r and 
from the CLV angle distribution [Fig.[26l^a,b,e)]. For this 
property, we call these modes tangled spurious modes. 
The subspace composed of them is also hyperbolically 
isolated from its complementary [red dashed line and 
squares in Fig. [TSTcd)]. One should be reminded that 
the difference between the tangled spurious modes and 
the physical ones lies in the sign of their Lyapunov ex- 
ponents; while the exponents for the physical modes can 
take any sign, those for the spurious modes are always 
negative, and thus they decay without intervening in the 
dynamics even if they are internally connected. The ef- 
fective dimension of the system is therefore always given 
by the number of the physical modes, A^ph, without con- 
tribution from the tangled spurious modes. 

The appearance of the tangled spurious modes is in 
fact, again, due to the discrete Laplacian operator of the 
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coupling term in Eq. (p!5]) . but it takes place in a different 
manner as the tent coupled-map lattice because of the 
continuous time. Here, if we take into account only the 
coupling term, its growth rate is simply given by the 
real part of its eigenvalues, specifically A = —{2/K^)[l — 
cos(27r/c/L)] for the eigenmodes e^^^^*/^. It is monotonic 
with respect to |/c|, since \k\ varies only up to I//2, and 
hence so is the peak wave number of the CLV power 
spectra in the stepwise region of the Lyapunov spectrum 
(Fig. [25l to be compared with Fig. [22]). The values of the 
Lyapunov exponent in this region are well approximated 
by A given above with the corresponding k (black line 
in Fig. [24]) . By contrast, for larger k or j, the gap in A 
between different k shrinks, as opposed to the case of the 
continuous Laplacian operator. This facilitates the order 
exchange of the Lyapunov exponents, i.e., the violation 
of the DOS, due to their intrinsic fiuctuations from the 
dynamics, and thus produces the tangled spurious modes. 
With increasing /i, both physical and tangled spurious 
modes increase their numbers, replacing the steps of the 
spurious modes in between (cf.. Fig. [M|) . This finally 
gives rise to the coalescence of the two groups; as soon as 
the last physical mode and the first tangled spurious one 
encounter, all the modes are connected and thus become 
physical modes. The case for h = 0.8 is shown in Fig. [271 
All the CLV angle distributions have finite densities at 
^ = and TT except for those concerning one or both of 
the two zero Lyapunov modes [Fig. [27r a)]. The overall 
connection of the Lyapunov modes is also confirmed form 
the violation of the DOS [Fig. [STfb)]. That is to say, 
with increasing /i, the number of the physical modes, 
or the inert ial manifold dimension as we conjecture, is 
suddenly doubled when the tangled spurious modes join 
the physical manifold. Although such a large value of 
h (around 0.7 in this case) corresponds to a somewhat 
unphysical situation when we view the CGL lattice as an 
approximation to the original CGL equation, the physical 
implication and consequence of this "crisis" remain to be 
clarified. 
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FIG. 28. (color online). Lyapunov spectrum for the 2D KS 
equation. Red dashed line shows the values of A*^^-* estimated 
from the dominating wave number (kx,ky) assigned trivially 
by the index (see text). Left inset: close-up of the begin- 
ning of the spectrum. Right inset: Lyapunov exponent A*^^-* 
vs /c^ — /c^, where k^ = k^ -\- ky is the square of the peak 
wave number for the power spectrum of the jth CLV (see 
Fig. [29)) . Blue crosses and black circles correspond to the 
physical and spurious modes, respectively. Red dashed line 
indicates A*^'^-* = k^ — k^. The Kaplan- Yorke dimension in this 
case is about 25.1. 
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V. 2D KURAMOTO-SIVASHINSKY EQUATION 
A. Definition and numerical scheme 

In the previous sections we have investigated the 
physical-spurious decoupling in ID dissipative systems. 
Given the generality of our finding, and also the general 
existence of a finite-dimensional inertial manifold in dis- 
sipative systems, one expects that the same decoupling 
should also take place in systems of higher dimensions. 
This is indeed the case, as we shall demonstrate in this 
section for the 2D KS equation as a generic example. 

The 2D KS equation is defined for a real-valued field 
u{x^ y^t) as [40| 



du 
'dt 



1 



V'u-{V')'u-^{Vu)', 



(16) 



FIG. 29. (color online). Power spectrum S{kx,ky) of the 
trajectory (a) and of the CLVs (b-f) for the 2D KS equation. 
The spectra are shown for kx,ky G [— /ccut, /ccut] with kcut = 
10 • 271 /L. 



with V = {d/dx,d/dy), x e [0,L^], and y e [0,Ly]. 
Here we adopt PBG u{x^y^t) = u{x-\-Lx^y,t) = u{x,y-\- 
Ly^t) with L = Lx = Ly = 24 for the sake of simplic- 
ity. For integration we use the same operator-splitting 
method as for the ID KS equation and the pseudospec- 
tral method with kcut = 10 • 27r/L, i.e., with 11 x 11 
Fourier modes. Integration is performed over a period 
of roughly 2 x 10^ after discarding a transient period of 
1200. 
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B. Results 

For L = 24 the 2D KS equation shows spatiotemporal 
chaos as in the ID counterpart studied in Sec. |TT1 This 
leads to a similar structure of the smooth region of the 
Lyapunov spectrum [left inset of Fig. [28l cf. Fig. [Hb)], 
which consists of physical Lyapunov modes carrying in- 
formation on the dynamics. As for the spurious region, 
despite its apparently intricate structure in contrast with 
the ID spectrum [see Fig. [TJa)], we find that the modes 
therein possess essentially the same trivial properties, as 
shown below. 

First, unlike the physical modes which show relatively 
broad power spectra similar to the dynamics [Fig. [SoTa- 
c)], each spurious mode has a spectrum S{kx^ ky) sharply 

peaked around a certain wave number k = \ k^ + k^ 



[Fig. [29r d-f)]. which is large and increases with the index. 
The value of its Lyapunov exponent A^-^^ is then essen- 
tially determined by the linear terms of the 2D KS equa- 
tion ([T6]) . specifically A^-^^ = /c^ — /c^, as indeed confirmed 
in the right inset of Fig. [28] (black circles). Moreover, as 
in the ID KS equation, the peak wave numbers (kx^ky) 
of the spurious modes are trivially determined by the ge- 
ometry: for a finite box of size L^xLy, the wave numbers 
can take values kx = 27rnx/Lx and ky = 27rny/Ly with 
intergers rix and Uy. Their values are then assigned to 
each mode in such a way that the resulting Lyapunov ex- 
ponent A^-^^ = k'^ — k^ decreases monotonically. Taking 
into account the multiplicity for 2D PBC, this trivial as- 
signment yields a spectrum indicated by the red dashed 
line in Fig. [28l which reproduces well the structure of the 
spurious region. 

To locate exactly the threshold between the physical 
and spurious modes, one needs to measure the hyperbol- 
icity properties as in all the other cases. Figure ISOT a) 
shows the time fraction Ur^''^^ of the DOS violation for 
arbitrary pairs (ji, J2)- Although internal connections 
are found for some neighboring steps of spurious modes, 
which are probably caused by mixing of several Fourier 
modes due to accidentally close values of the Lyapunov 
exponent [see Fig. [29r d-f)]. an unambiguous threshold is 
found at A^ph = 121, for which ui-^''^^^ = for ah the 
pairs with ji{J2) < A^ph < J2{ji) even at the smallest r 
we used [white dashed line in Fig.lSOTa)]. Varying r does 
not unravel the internal connections of spurious modes 
presumably, as u^'^'^ decays in the fastest case expo- 
nentially [Fig.lSOTb)]. whereas it would decay faster for a 
hyperbolically isolated pair [see Figs. flST a) and flHT a)]. 

The same conclusion is reached from the angle distribu- 
tions pv((9) of the CLVs [Fig.[30i;c,d)]. Pairs of CLVs for 
physical modes can take any angles and experience rather 
frequent tangencies, while the angle distribution for spu- 
rious pairs is peaked around 7r/2 and decays quickly in 
the tail. A closer look at this tail reveals that hyperbol- 
ically isolated pairs, or those which fulfill the DOS, de- 
cay as pv{0) ~ exp(— const. /^) or faster [two innermost 
plots in Fig. ISOTd)] in contrast with the slower decays 
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FIG. 30. (color online). Hyperbolicity properties for the 2D 
KS equation, (a) Time fraction z/^ ^''^^ of the DOS violation 
for arbitrary pairs with r = 0.2. The black color indicates 
^ui>J2j _ Q^ j^^ hyperbolically isolated pairs. The white 
dashed lines indicate the threshold Aph = 121 between the 
physical and spurious modes, (b) Time fraction Ur^ of the 
DOS violation for pairs of neighboring Lyapunov exponents as 
a function of r. (c) Angle distributions pv{0) for neighboring 
CLVs of indices (89,90), (129, 130), (109, 110), (121, 122), and 
(137, 138) from outermost to innermost curve. Notice that 
they are not in order with respect to the index, (d) Pv(0) 
against 1/0. The ordinate is averaged over both sides of the 
distribution. The order of the curves is the same as in (c) 
from upper right to lower left, (e) Angle distribution ps(0) 
between the physical subspace (1 < j < Aph) and the spurious 
one (Aph + 1 < i) with Aph = 121 (green solid line). The red 
dashed line shows the angle distribution ps((/>) between the 
subspaces of indices 1 < j < 109 and 110 < j. (f) The 
same plots as (e) but shown against 1/0 (green diamonds 
for the physical and spurious subspaces, and red squares for 
the other). Each color corresponds to the same index pair 
throughout in the panels (b-f). 



found for pairs of the physical modes (black circles and 
red squares) or the internally connected spurious modes 
(blue triangles). The hyperbolic decoupling at A'ph = 121 
is established by the same fast decay of the angle distri- 
bution Ps(</>) between the two subspaces defined thereby 
[green full curve and diamonds in Fig. [SOf e.f)]. Defin- 
ing the threshold A'ph at a false value leads to a clearly 
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slower decay in Ps(^) [red dashed curve and squares, ob- 
tained with the threshold corresponding to the red line 
and squares in Fig. ISOTb-d)]. In short, the hyperbolic 
decoupling between the physical manifold and the spu- 
rious modes takes place in the same way as for the ID 
systems, with the peak wave numbers for the spurious 
modes distributed here in 2D Fourier space. 



VI. DISCUSSION 

In the present paper, we have shown the hyperbolic 
decoupHng of the tangent space in spatially-extended dis- 
sipative systems into a finite set of A/ph physical modes 
characterizing all the physically relevant dynamics, and 
a remaining set of spurious modes, which represent triv- 
ial, exponentially decaying perturbations. Physical and 
spurious modes are characterized by the existence and 
absence, respectively, of tangencies between the associ- 
ated CLVs. In contrast to physical modes which are 
densely connected to each other through frequent tan- 
gencies, the spurious modes are hyperbolically isolated 
from all the physical modes, and therefore solely decay 
to zero without activating any physical modes. The exis- 
tence or absence of tangencies between Lyapunov modes 
can be investigated both from the CLV angle distribution 
and the violation of DOS. One can also take into account 
linear combinations of physical or spurious modes by 
measuring the angle between the two subspaces spanned 
by the physical and spurious modes, respectively. This 
leads to the unambiguous determination of the number of 
the physical modes, A^ph- We conjecture that this num- 
ber corresponds to the minimal dimension of the iner- 
tial manifold. Moreover, we have shown that the hy- 
perbolic decoupling between the physical and spurious 
modes takes place quite generically in dissipative sys- 
tems, as we have found evidence in PDEs in one or higher 
spatial dimensions, as well as lattice systems with con- 
tinuous or discrete time. 

A few remarks about related earlier studies are in or- 
der. As already mentioned, the stepwise structure of the 
Lyapunov spectrum in the spurious region was already 
reported by Garnier and Wojcik for CGL amplitude tur- 
bulence, leading them to the conjecture that these are 
inactive Lyapunov exponents [36|. This is pertinent as 
we have revealed in the present study, but one should be 
reminded that the stepwise Lyapunov spectrum does not 
define or indicate the spurious modes. The steps do not 
form at a single well-defined threshold index: they are 
actually slightly inclined, and get fiatter and flatter as 
one goes deeper into the spurious region. As a matter 
of fact, the physical region typically extends to the first 
few steps of the spectrum: for the ID KS system studied 
here the steps begin at j = 40 but the spurious modes 
are from j = 44. This gap can be quite large; it is found 
to be 14 for the coupled tent maps studied in Sec. IIV A[ 
The formation of the steps is only a consequence of the 
sharp power spectrum and the choice of the boundary 



conditions. 

In a different work, Ikeda and Matsumoto studied the 
projection of the dynamics onto Fourier modes or on 
the Gram-Schmidt vectors associated with the Lyapunov 
modes [41]. They measured in particular information- 
theoretic quantities and argued that information does 
not flow from a high wave number region to a lower 
wave number region. However, similarly to Lyapunov ex- 
ponents. Gram- Schmidt vectors or Fourier modes alone 
cannot characterize the separation between physical and 
spurious modes, and thus cannot provide the exact in- 
dex of the threshold. It is the CLVs that bear the true 
physical properties of tangent space and thus are able 
to define unambiguously the physical and spurious Lya- 
punov modes from their hyperboHcity properties. 

Concerning the difference between the CLVs and the 
Gram- Schmidt vectors, while it is clear that properties 
directedly linked to hyperbolicity such as angle distribu- 
tions can only be captured by CLVs, one may formally 
study DOS using either CLVs or Gram-Schmidt vectors. 
However, clearly, the finite-time Lyapunov exponents de- 
fined using the two sets of vectors are a priori different 
for any finite times. Indeed, although these two defi- 
nitions seem to give qualitatively similar results [42], it 
remains true that the mathematical connection to hy- 
perboHcity is only proved for CLVs [30|, ^M- Thus, in 
the absence of a detailed study on this point, one cannot 
avoid computing CLVs for drawing any firm conclusion 
on hyperbolicity properties of dynamical systems. 

We should also mention that all the angle distributions 
and the finite time Lyapunov exponents used for the de- 
termination of DOS are not invariant against coordinate 
transformations. The splitting, however, between phys- 
ical and spurious modes is coordinate independent be- 
cause by smooth transformations a nonzero angle cannot 
be made to vanish and vice versa. Furthermore, the type 
of the essential singularity found in the angle distribution 
of critical pairs is invariant against smooth transforma- 
tions of the coordinates, as can be seen by considering 
how the angles are transformed with the coordinates. A 
simple example is provided by the two versions of the KS 
equations, Eqs. ([T]) and ([IS]), which are connected by a 
linear transformation [40]. These considerations further 
emphasize the generality of our results. 

Pursuing our conjecture that physical modes form the 
local approximation of the inertial manifold, it is cru- 
cial to seek a connection between CLVs and phase-space 
properties, and in particular to find a way to construct 
the inertial manifold using information of the physical 
CLVs. A straightforward construction is a manifold lo- 
cally defined at each point of the trajectory by the phys- 
ical CLVs, but strictly speaking, one also has to consider 
the nonlinear growth of finite-size perturbations, which 
is obviously not captured by Lyapunov analysis. Because 
of this, it is still in principle possible to consider a system 
in which hyperbolically isolated modes correspond to ac- 
tive degrees of freedom, such as a lattice of uncoupled 
hyperbolic maps as a trivial instance. However, given 
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that the generic systems we studied are by no means hy- 
perbohc, it is unlikely that that type of decoupling has 
any relation to the hyperbolic decoupling studied in the 
present paper. 

Our results are also important from a practical point of 
view. Since the physical properties of the PDEs at stake 
are carried by the physical modes, a faithful numerical in- 
tegration needs to incorporate at least as many degrees of 
freedom as the physical modes. In other words, the num- 
ber of the physical modes, A^ph, sets a lower bound for 
the spatial resolution of the simulations. Moreover, the 
extensivity of A^ph implies that measuring it in systems of 
moderate size suffices to have reliable estimates for arbi- 
trarily large system sizes. Our results may also be useful 
for controlling spatiotemporal chaos, in particular, when 
determining the minimal number of constraints for full 
control of such a system. A possible application in this 
context would be the suppression of ventricular fibrilla- 
tion, the major reason of sudden cardiac death, for which 
promising schemes have been proposed recently applying 
low amplitude electric shocks over grids of points or lines 
[43( . The number of the physical modes would then help 
determine an appropriate grid resolution. Note however 
that, having no access to the basis constituting the in- 
ertial manifold, one may actually need more degrees of 
freedom or constraints than A^ph to achieve faithful simu- 
lations or full control of a dissipative system. It is there- 
fore of both fundamental and practical importance to 
elucidate the structure of the inertial manifold in phase 
space in the fine of the present study. 

In summary, we have shown, using Lyapunov anal- 



ysis and, in particular, CLVs, that the tangent space 
of spatially-extended dissipative systems can be divided 
into two parts: a finite-dimensional manifold spanned by 
the physical modes bearing all the physical information 
and the remaining set of spurious, exponentially decaying 
modes. We have demonstrated that the spurious modes 
are hyperbolically separated from all the physical modes 
and hence satisfy DOS. The separation holds for any lin- 
ear combination of the physical or spurious modes, as the 
minimal angle between the physical and spurious sub- 
spaces is bounded away from zero. The number of physi- 
cal modes has been interpreted as the number of effective 
degrees of freedom needed to describe the dynamics. This 
leads to the conjecture that the physical modes constitute 
a local approximation of the inertial manifold. We hope 
that these results will trigger mathematical studies in or- 
der to put these issues on a more rigorous basis. It is also 
interesting to perform further numerical investigations of 
other dissipative systems such as the Navier- Stokes equa- 
tion, for which no rigorous proof of the existence of an 
inertial manifold is known. 
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